function H = generate_channel(Nt, Nr, L)  %(NLOS)

AOD = pi*rand(L, 2) - pi/2;  %-2/pi~2/pi AOD代表发送的方位角和仰角 相当于生成了L×2 个-0.5 π ~0.5 π的角
AOA = pi*rand(L, 2) - pi/2;  %-2/pi~2/pi
alpha = (randn(1,L)+1i*randn(1,L))/sqrt(2);%信道的衰落系数
%为了区别LOS径和NLOS径，我们默认第一径是LOS径，给了归一化的能量为1。
% 而NLOS径的能量则是高斯变量， 并且为了体现出与LOS径的能量差， 给了一个衰减系数。
H = zeros(Nr, Nt);
for l=1:1:L
    ar = array_response(AOA(l,1),AOA(l,2), Nr);
    at = array_response(AOD(l,1),AOD(l,2), Nt);
    H = H + sqrt(Nr * Nt/L)*alpha(l)*ar*at';
end